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We investigate chaos in mixed-phase-space Hamiltonian systems using time series of the finite¬ 
time Lyapunov exponents. The methodology we propose uses the number of Lyapunov exponents 
close to zero to define regimes of ordered (stickiness), semi-ordered (or semi-chaotic), and strongly 
chaotic motion. The dynamics is then investigated looking at the consecutive time spent in each 
regime, the transition between different regimes, and the regions in the phase-space associated to 
them. Applying our methodology to a chain of coupled standard maps we obtain: (i) that it allows 
for an improved numerical characterization of stickiness in high-dimensional Hamiltonian systems, 
when compared to the previous analyses based on the distribution of recurrence times; (ii) that 
the transition probabilities between different regimes are determined by the phase-space volume 
associated to the corresponding regions; (iii) the dependence of the Lyapunov exponents with the 
coupling strength. 
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I. INTRODUCTION 

In weakly chaotic Hamiltonian systems regions of regu¬ 
lar (periodic and quasi-periodic) and chaotic motion typ¬ 
ically coexist in the phase-space Hi- In high dimen¬ 
sions, due to Arnold diffusion, all initial conditions lead¬ 
ing to chaotic motion are connected in the phase-space 
building a single chaotic component [l|. Even if the vol¬ 
ume of the regular regions becomes vanishingly small, as 
expected for high-dimensional nonlinear systems, the dy¬ 
namics inside the chaotic component of the phase-space 
is strongly affected by such regions. This happens be¬ 
cause trajectories approaching non-hyperbolic regions or 
regular motion remain a long time close to them before 
visiting again other parts of the chaotic component of 
the phase-space. This signature of weak mixing (or weak 
chaos) is known as stickiness m. 

Since Chirikov-Shepelyansky Q, the main quantifi¬ 
cation of stickiness in Hamiltonian systems has been 
through the fat-tail distribution of Poincare recurrence 
times (see, e.g., 0, i i). An alternative approach 

is to use finite-time Lyapunov exponents (ETLEs) [13- 
[T3| . with recent applications using lar ge d eviation tech¬ 
niques 0, [l3| and the cumulants [13,113] of the ETLE 
distribution. In area-preserving maps, stickiness generi- 
cally occurs at the border of 2-dimensional Kolmogorov- 
Arnold-Moser (KAM) island [I[ {i.e., at 1-dimensional 
tori). The recurrence time is a measure of the time the 
trajectory spends around such structures before return¬ 
ing to the chaotic sea (stickiness happens also to one- 
parameter families of parabolic orbits EMa and even 
to isolated parabolic fixed points Near the 

non-hyperbolic structures, the local instability of chaotic 
trajectories is reduced so that ETLEs can be used to char¬ 
acterize phase-space regions of interest [iMiiio, mi. 


Stickiness has been studied also in higher-dimensional 
systems 0, [Hill El, mi. Ell : recurrence times can 

be due to different non-hyperbolic regions and tori of dif¬ 
ferent dimensionalities [23|. An improved characteriza¬ 
tion of stickiness events (long recurrence time) requires 
thus to measure the number of stable and unstable di¬ 
rections in the trajectory during this event. Eroeschle 
conjectured that lower-dimensional tori could not ex¬ 
ist El El El- In early studies in the 80’s such events 
of stickiness to lower dimensional tori were reported in 
some systems [^1 but were not found in other exam¬ 
ples Eoj. Even if invariant tori do not exist, small lo¬ 
cal Lyapunov exponent could effectively act as a lower¬ 
dimensional trap. This is similar to almost invariant sets 
[H, El , which are regions in phase-space where typical 
trajectories stay (on average) for long periods of time. 

In this paper we introduce a methodology that uses 
time-series of local Lyapunov exponents to define regimes 
of ordered, semi-ordered and totally chaotic motion and 
obtain an improved characterization of stickiness in high¬ 
dimensional Hamiltonian systems. We illustrate this gen¬ 
eral procedure in a chain of coupled standard maps and 
confirm that stickiness events of different times length 
are dominated by trajectories with different ETLEs. A 
significant improvement of the characterization of sticky 
motion in high-dimensional systems is found. We also 
characterize the ETLEs for small couplings and compare 
them to expected universal properties in fully chaotic sys¬ 
tems [llj. The method proposed here is general and can 
be used to investigate Hamiltonian systems in any di¬ 
mension. 

The paper is divided as follows. In Sec. mi we describe 
the Hamiltonian model we use to illustrate our method. 
In Sec. m we introduce our method to compute and 
analyze time series of local Lyapunov exponents. This 
methodology is then applied to the symplectic model of 
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coupled standard maps in Sec. lIVI Section IVl summarizes 
the main results of the paper. 

II. THE COUPLED MAPS MODEL 

We use a time-discrete 27V-dimensional Hamiltonian 
system obtained as the composition T o M of indepen¬ 
dent one-step iteration of N symplectic 2-dimensional 
maps M = (Ml,..., Mat) and a symplectic coupling 
T = (Ti,..., Tat). As a representative example of 2- 
dimensional maps we choose for our numerical investi¬ 
gation the standard map: 

M. I I _ I Pi + KiSin{2Trxi) mod 1 \ 

\ Xi j \ Xi + Pi + KiSm{ 2 Trxi) mod 1 j ’ 

and for the coupling 



with ^ij = ^j^i = (all-to-all coupling). The moti¬ 

vation for working with this system is that in the limit 
of small coupling ^ 0 it can be understood looking 

at the dynamics of the N uncoupled maps. This system 
was studied in Refs. [1,[3 using recurrence time distribu¬ 
tion. This allow us to critically compare the benefits of 
our methodology. In all numerical simulations we used 
Ki = 0.5214 for the map Mi and K 2 = K 3 = 0.5108 for 
the maps M 2 and M 3 . 

III. METHOD 

In this section we describe the method proposed in 
this work. To be illustrative, we present numerical sim¬ 
ulations for the system defined in Sec. [Ill 

A. Lyapunov spectrum and the classification of 
ordered, semi-ordered or semi-chaotic, and chaotic 
regimes 

Consider a chaotic trajectory in a closed Hamiltonian 
system which, after reducing the phase-space dimension 
due to global invariant of motion, has N degree of free¬ 
doms. For long times t the trajectory ergodically fills 
the whole chaotic component of the phase-space which 
is characterized by a spectrum of N Lyapunov expo¬ 
nents jy}, where > A^'^^ ..., A^^ > 0 [t^ . 

The central ingredient of our analysis is the spectrum 
of FTLEs computed along a trajectory during a win¬ 
dow of size u! where we obtain a time dependent spec¬ 
trum {A-“^}(t) = {A-“^}. The window size w has to be 
sufficiently small to guarantee a good resolution of the 


temporal variation of the A^^^^’s, but sufficiently large in 
order to have a reliable estimation (see Refs. [HTI - Ii^ L 
The probabilit y d ensity function of a|"^^ has been exten¬ 
sively studied [3, Hol - H^ . Here we go beyond the study 
of the probability density function and explore temporal 
properties in the time series of {a[“^}. 

Figures dKa) and (b), for iV = 2 and 3 respectively, 
show the time series of a[“\ (i = 1,..., A^). The sharp 
transitions towards « 0 motivates the classification 
in regimes of motion | 2 l| as (a) ordered (A ^‘(^2 ~ 0 ); 
(b) semi-ordered or semi-chaotic (A^*^^ > 0 ; A^*^^ ~ 0 ); 
and (c) chaotic (A ^‘)^2 > 0). For a system with N de¬ 
grees of freedom we will say that the trajectory is in a 
regime of type if it has M local Lyapunov exponents 
A^”^^ > Ei, where £i <C A-°“^ are the small thresholds. 
This means that and are ordered and chaotic 
regimes respectively. Whenever there is no ambiguity, we 
will drop the superscript = Sm to have a simpler 

notation. 

Practical implementations of the general method de¬ 
scribed above require the choice of a few parameters and 
conventions. First of all, the window size w and the 
threshold £i directly affect the classification in regimes. 
They can be thought as the phase-space resolution of the 
analysis and should be chosen so that it provides max¬ 
imal information about the regions of interest. Unless 
stated otherwise, we use lu = 100 and Ei « 0.10(A-‘^^), 
where (...) denotes average over t, where t = 1 ,..., ti 
(even though the classification in regimes is strongly oj- 
dependent, our conclusions are not sensitively affected 
by variations around the chosen values). Another impor¬ 
tant choice is the method for computation of the FTLEs. 
We use Benettin’s algorithm [^, , which includes the 

Gram-Schmidt re-orthonormalization procedure. The 
decreasing order of a[“^ is valid on average, but inver¬ 
sions of the order (A^!))^]^ > A,-*^^) may happen for some 

times t and we have chosen to impose the order of A^^^^ 
for all t. Finally, it is possible to decide how to sample 
the time series A^'^^ While the FTLEs are defined for 
all t, there is a trivial correlation between the values of 
FTLEs inside a window of size to because they are com¬ 
puted using the same points of the trajectory. In order 
to avoid this trivial correlation the series of A^"^^ can be 
computed using non-overlapping windows, i.e. plotting 
A^'^^ only every w time steps (a choice we adopt in our 
simulations). 


B. Identifying phase-space regions 

In order to understand the properties of the time se¬ 
ries {A-“^} it is useful to consider the phase-space regions 
associated to each regime Sm- We denote by p[A) the 
phase-space volume (Liouville measure) of region A in the 
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Figure 1. (Color online) Illustration of the method proposed to define the regimes Sm- Time series of the spectrum of FTLEs 
with i = 1,..., A^, for the map ([l])-© with ^ = 10~®. Panel (a): Case N = 2 and the thresholds ei = 0.1 and 
£2 = 0.05 are represented by dash-dotted and dotted lines respectively. Panel (b): Case N = S and the thresholds £i = 0.1, 
£2 = 0.08 and £3 = 0.06 are represented by dash-dotted, dotted and dashed lines respectively. 


bounded phase-space F, i.e. ^i{T) = 1. The most impor¬ 
tant distinction is between the regions of regular Fj-eguiar 
and chaotic Fchaos motion. In Hamiltonian systems, typ¬ 
ically /i(rchaos) > 0 and ir(ri.eguiar) > 0. In principle, 
the regular region Fj-eguiar can be subdivided according 
to the dimensionality of the tori. However, according 
to Froeschle’s conjecture, in a 27V-dimensional phase- 
space, tori with dimension N have positive measure and 
thus M(rreguiar) = ^(rtori) [1 [H, Si] . For TV > 1, the 
chaotic region Fchaos is expected to build a single ergodic 
component because tori of TV dimension do not parti¬ 
tion the 2TV-diinensional phase-space in different regions 
and therefore any chaotic trajectories eventually explores 
(through Arnold Diffusion) the whole Fchaos- Our inter¬ 
est is not to test the Froeschle conjecture or Arnold diffu¬ 
sion, but to show the insights about the chaotic dynam¬ 
ics we can obtain using the time series of together 

with the definition of the regimes Sm- One application is 
to use the regimes Sm to split the chaotic component of 
the phase-space in meaningful components. This is done 
by considering the set of points in the phase-space 

leading to each regime Sm as 

= lim xtixt G Sm), (3) 

tL—^OO 

where is the total length of the trajectory and Xt G Sm 
indicates that at time t the trajectory at Xt had G 

Sm- 

Figure [5] shows numerical estimates of the phase-space 
regions obtained for each regime Sm in the chain of cou¬ 
pled maps defined in Sec. m The regime Sq (or the or¬ 
dered regime) is associated to region localized close to the 
border of the KAM island of the uncoupled case (com¬ 
pare to Fig. [2^). Points which belong to the regime Si are 
closer to the center of the torus from the uncoupled case. 
This suggests that when trajectories are inside the region 
related to regime , they more likely penetrate inside the 


torus from the uncoupled case. In the chaotic sea both 
regimes Si and S 2 are visible. These results are natu¬ 
rally understood in the perturbative limit (small coupling 
^ <C 1). The regime Sq corresponds to A-^™^ « 0 for ev¬ 
ery T = 1 ,..., X, which is expected when the trajectory 
is stuck close to the X-dimensional tori built as the prod¬ 
uct of the 1-dimensional tori of the uncoupled maps. In 
contrast, Sm for TVf > 0 implies that at least one FTLE 
A-”^^ 0 and therefore the trajectory projected in one 

map can be both in the chaotic and regular regions (e.g.. 
Si for X = 2 can be obtained from Aj“^ » 0, ~ 0 

or from A^"^^ ~ 0, X^'’ :$> 0). Altogether, these obser¬ 
vations confirm that our method allows for a meaningful 
division of the chaotic component of the phase-space and 
can thus be used to identify regions of interesting dynam¬ 
ics. In the case partial barriers exist inside the chaotic 
component - such as in area-preserving maps with mixed 
phase space [^ - we expect the regions obtained through 
our method to depend weakly on uj and to coincide with 
those obtained from the partial barriers. 


IV. RESULTS 

In this Section we apply the Lyapunov time-series 
methodology described in Sec. lHII to the 2X-dimensional 
system defined in Sec. [Ill We compute and interpret four 
basic properties of the method: the total time spent 
in each regime (residence time), the transition between 
regimes, the consecutive time in each regime, and the 
scaling of Lyapunov exponents. 


A. Residence time in each regime 

The first and most basic quantity we measure is the 
probability P{Sm) of finding the trajectory in each 
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Figure 2. (Color online) Phase-space projected in {xi,pi) for different configuration of the A'^-coupled standard maps defined in 
Sec. HD (a) = 1 (uncoupled case), showing 10 ^ randomly started initial conditions and plotting as dots 10 "^ iterations of each 

of them. A large KAM island can be seen at the center of the plot; (b) N = 2 and coupling strength ^ = 10“®; (c) N = 3 and 
f = 10“®. Symbols with different colors in (b,c) show points xt £ Sm belonging to regimes So (blue circles), Si (red points), 
and S 2 (green points). These points were computed starting a single trajectory in the chaotic region of all maps and iterating 
it 5 X 10® times. 


regime, defined as the fraction of the total time tp 
that Xt e Sm (*-e. P{Sm) = ^teSM where 
SteSM = 1 if t G Sm and SteSM = 0 otherwise). 



Figure 3. (Color online) Residence time in each regime Sm- 
(a.) N — 2 with ei = 0.1 and £2 = 0.05. (b) N = 3 with 
£1 = 0.1, £2 = 0.08, and £3 = 0.06. In (a) the values obtained 
with to = 100 are compared (gray curves) with results for 
a) = 50 and lj — 500. Only for the case M = 1 the gray 
curves (right for a; = 50 and left for uj = 500) show a shift in 
the x-axis Estimations for each ^ are based on a trajectory 
with length ~ 10 ^®. 


Figure [3] shows the probabilities P{Sm) for the map 
with = 2 ,3 as a function of the coupling strength 
We now explain the behavior of P{Sm) with ^ by 
discussing the effect of coupling ^ on the phase-space re¬ 
gions associated to Sm, as defined in Sec. IIIIBI By the 
ergodicity of Fchaos, P{Sm) corresponds to the (normal¬ 
ized) volume of the region related to regime Sm in the 
phase-space 


P{Sm) 


^J■{SM) 

M(f^chaos) 


KSm) 

1 - ^(Ftori) ’ 


(4) 


The results of Fig. [3] show that the chaotic region is the 
largest region in phase-space for any coupling, while the 
region associated to Si has a larger volume than Sq for 
couplings ^ < 1.3x 10“®. For larger ^ we see oscillations 
with a local maximum close to ^ 2 x 10 “® for the cases 

M = 0 and M = 1 . 

We now interpret the ^ dependence observed in Fig. [3] 
by arguing how the different terms in Eq. o vary with 
We denote by /i(L'j) the measure of tori for the j-th 
map with control parameter Kj in the uncoupled case 
^ = 0 (which we assume to be approximately equal to 
the measure of the KAM islands). For small coupling 
^ « 0 we expect that most tori of the uncoupled maps to 
survive and therefore: 

• ^(rtori) ~ n^i which in the simple case of 

l^{Uj) = pL for all j reduces to /i(Ftori) ~ ■ 


• p,{Sm=o) corresponds to a small volume around 

r tori 5 e. fJ-{SM=o) ^(rtori) ~ P{^o)/{^ + 

P{So))- 


• For p,{Sm^o) we have that N — M maps are in 
their corresponding KAM island (with probability 
p{Uj)) and M maps in the chaotic area (with prob¬ 
ability 1 — /r(C/fc)). For example, for iV = 3 and 
M = 2 we have that 


fJ.{SM=2) = /i(Cfi)[l - fJ-{U2)][l - fJ-iUs)] 



























+ ^^{U2)[i - KUiW - KUs)] 

+ ^{Us)[l - ^iiU^)][l - fi{U2)]. 
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In general this leads to 

KSm) n n 

where the last product is over all j 

except j € {ji, ■ ■ ■ jJm}- In the simple case of 

/r(t/j) = fj,, it reduces to ^i{Sm) ~ ~ 

We now consider the effect of growing In the spirit 
of the KAM theorem, the tori of the coupled maps (gen¬ 
erated as the product of the N maps) are expected to be 
robust to small couplings which act as a perturbation. 
This explains why the curves in Fig. [3] are essentially 
flat for small Increasing ^ even further, the nonlin¬ 
earity of the system increases and therefore /r(rtori) is 
expected to decrease (/r(rtori) —t 0 for ^ ^ 0). This 
reduction of the tori leads to an increase in the denom¬ 
inator of Eq. ([3]) and explains the observed tendency of 
reduction of P{Sm) for all regions related to stickiness 
(M < N). Indeed, for > 0.5 no signature of tori or 
stickiness was detected numerically and P{Sm=n) = 1 - 
The nontrivial dependencies of P{Sm<n) in Fig. [3] ap¬ 
pear at^~2xl0 “^ values, close to the values of ^ for 
which the last tori disappear (see also Fig. 7.2 in Ref. @). 
In this regime the volume of the tori is already negligible 
M(rtori) ^ 0 but stickiness is still effective (notice that 
even zero measure non-hyperbolic sets can lead to stick¬ 
iness [3, [13 )■ The denominator in Eq. dl is therefore 
1 — ^(Ftori) 1, not significantly affected by further in¬ 
creases of and therefore not driving the reduction of 
PiSM<N)- Small variation of a control parameter of the 
system (in this case are known to lead to sensitive 
creation and destruction of tori, with non-trivial depen¬ 
dencies on the stickiness [3- We can thus expect that - 
close to the disappearance of the tori - the small volume 
of stickiness regions h{Sm<n) to fluctuate with ^ lead¬ 
ing even to an increase with It is interesting to note 
that this non-trivial increase with ^ appears for P{Sm=o) 
in Fig. [3] precisely when the curves P{So<m<n) show a 
sharp decreasing fluctuation. This suggests an exchange 
between measure of different sticky regions associated to 
regimes Sm<n^ without interference of the much larger 
fully chaotic component Sm=n- 

B. Transitions between regimes 

We now focus on the transition between regimes. The 
simplest analysis correspond to the two-time (joint) prob¬ 
ability P{Sm Sm’), computed as the fraction of the 
total trajectory time that Xt G Sm and Xt+i G S'm'- 
The probabilities considered in the previous section can 
be obtained as J2 sm t Sm') = P{Sm') and 

P{Sm —t Sm') = P{Sm)- Figure Ua) shows the 



Figure 4. (Color online) Transition between regimes Sm as a 
function of the coupling strength (a) Transition probability 
P{Sm —>■ Sm')', (b-d) Conditional probability Pm,m' defined 
in Eq. © of moving to M' given that the trajectory was at 
M. Estimations for each ^ based on a single trajectory with 
length tL = 10^° in the case of A = 2 coupled maps and 
£i = 0.1 and £2 = 0.05. 

dependence of P{Sm —t Sm') on ^ for our model. We 
notice that P{Sm —t Sm') is equal to P{Sm' —t Sm)- 
This is expected considering that the system is ergodic, 
volume preserving, and time-reversible. The dependence 
of P{Sm Sm>) on ^ follows a similar pattern observed 
for P{Sm) in Fig. [31 More information is obtained from 
the conditional probability 

Pm,m' = P{Sm —t Sm'\Sm) = — ^ (5) 

which quantifies the probability that trajectories at Sm 
will move to Sm' ■ The results shown in Fig. |3](b-d) show 
for all Sm that (i) persistence in the same Sm {P[Sm 
Sm\Sm)) is dominant and (ii) the most likely transitions 
occur between neighboring regimes (e.^., P 2 ,i > P 2 ,o)- 
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The only (slight) deviations of this picture happen for 
large values of close to the disappearance of the KAM 
island. Altogether, these results confirm that in the per¬ 
turbative regime <C 1) stickiness happens approaching 
the region of regular motion of different maps one after 
the other (in opposite to a direct approach from So to 
Sm=n)- 


C. Consecutive time in each regime 


The results of the previous section confirm that resi¬ 
dence in the same regime is the dominant behavior. This 
motivates us to study the time tm spent consecutively in 
a regime Sm (*.e., tm is the time between two consec¬ 
utive transitions between different regimes, the first to 
Sm and the second out of Sm)- In a trajectory of length 

we collect a series of tm- We are mainly interested in 
the probability distribution P{tm) (or, equivalently its 
cumulative Pcum = =tm different S'm’s 

in the limit ^ oo. These distributions should be com¬ 
pared to the distribution of recurrence times r, defined as 
the time between two successive entries to a pre-defined 
recurrence region (usually taken in the fully chaotic com¬ 
ponent of the phase-space). Events in the tails of P{t) 
are associated to times for which the trajectory is stuck 
to the non-hyperbolic components of the phase-space and 
P{t) is the traditional method to quantify stickiness in 
Hamiltonian systems d, 0, H, 0 . 

The numerical simulations in Fig. [5] confirm that the 
distribution obtained summing Pcum{TM) for ordered 
and semi-ordered regimes (or Sm<n) is equivalent to 
cumulative distribution Pcum{T) obtained using recur¬ 
rences. This is in agreement with the association of 
long consecutive times in regimes of ordered and semi- 
ordered motion to long recurrence times. Looking at 
the individual distributions Pcum{TM) provide valuable 
additional information on the sticky motion. For semi- 
ordered motion (when 0 < M < N) we observe an ex¬ 
ponential tail after an intermediate decay with scaling 
j3 ~ 0.5. This behavior confirms the interpretation given 
in Ref. Q. More interestingly, the M = 0 case shows 
an asymptotic algebraic decay which characterize stick¬ 
iness. While the scaling is compatible with the results 
obtained using recurrence time, T’ciim('rM=o) obtained in 
our methodology provides a better characterization of 
the scaling (over several orders of magnitudes) and al¬ 
lows for an independent analysis of the different regimes. 
These properties are essential when dealing with high¬ 
dimensional systems (which may contain different pre- 
asymptotic regimes) and for an accurate estimation of 
the stickiness exponent 7 . Finally, panel (a) in Fig. O 
shows that all decays discussed above remain (qualita¬ 
tively) the same for different choices of lo, with the curve 
for M = 1 showing the largest sensitivity on w (as in 
Fig.IIa). 




r r 


Figure 5. (Color online) Comparison between our method 
and the analysis based on recurrence time. The cumulative 
distribution PcumiTM) of times tm is shown for each regime 
Sm<n for uj = 100 and (a) A = 2 and (b) A = 3. In 
(a) the gray curves show results for a; = 50 and uj = 500. 
Only for the case M = 1 the gray curves (left for oj = 50 
and right for oj = 500) show a shift in the a:-axis The 
cumulative distribution Pcum.{T) for recurrence times r to a 
region in the chaotic component of the phase-space (in Sm^n) 
for (c) A = 2 and (d) A = 3. For comparison, in panels 
(c) we show the results obtained combining the normalized 
curves for Mo + Mi (blue dotted line: divided by 1.7 x 10® for 
convenience of scale) of panel (a), and in (d) the normalized 
curves for Mo + Mi -|- M 2 (blue dotted line: divided by 10 ®) 
of panel (b). Results obtained using maps of Sec. ITIl with 
5 = 10~®, £1 = 0.1 and £2 = 0.05 for the case A = 2 and 
£1 = 0.1, £2 = 0.08 and £3 = 0.06 for the case A = 3. 


D. Scaling of Lyapunov exponents 


So far we have focused at the temporal properties of 
the time series of FTLEs and how they change with 
the coupling strength We now consider how the value 
of the Lyapunovs respond to an external perturbation, 
which in our case is the coupling to the other maps. It 
is known that the largest exponent is extremely 

sensitive to perturbation. More specifically, Daido’s rela¬ 
tion [ 2 ^, [s^ states that for small couplings ^ to another 
chaotic system, a universal logarithmic singularity is ob¬ 
served, 


^(o;—J-cso) ^(o;—j-cso) 


C 


( 6 ) 


|ln(OI’ 
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where are the unperturbed Lyapunov exponents 

and c is a constant and i = 1,... ,iV. This relation is 
valid for totally chaotic systems and a small mismatch 
between Lyapunov exponents of the uncoupled systems 
compared to their fluctuations [32|. Here we investigate 
the relation as a function of for distinct 

values of w and different regimes Sm- To this end we 
compute the temporal averages of the FTLEs for 

times t such that S Sm- 



(b) 




uj = 10^ 

- 



- 

(Af - 






■ 0.00 0.05 0.10 0.15 0.00 0.05 0.10 0.15 

lln(OI'' lln©!-' 


Figure 6. (Color online) Sensitivity of the FTLEs for small 
couplings. The difference between the finite-time Lyapunov 
exponent in the coupled ((A^^^)) and uncoupled ((A^“^)) maps 
as a function of |ln(^)|“^, where ^ is the coupling strength. 
Results are shown for = 2 (i = 1 , 2) and different time 
windows uj (a)-(b) 10 ^, (c)-(d) 10 "^, and (e)-(f) 10 ®. Black 
dashed lines in (c)-(f) are the expected linear behaviour, con¬ 
sistent with Eq. Panels in the left column (ace) were 

computed for the full time series, while on the right column 
(bdf) only FTLEs in the regime S 2 were used. The different 
colors correspond to different choices of threshold imposed to 
define the FTLE: (Ai“^)ei uses ei = 0.1(A;“^) while 
uses £2 = 0.9(Ai“^), where (A^^^) is computed over the full 
time series (left column). 


Our numerical simulations reported in Fig. [5] show 
that small values of ui lead to a situation in which 
ss at a finite value of ^ (Figs. EKb) and 

(d)), while larger values of w lead to situations in which 

C- These results depend cru¬ 
cially on our choice to impose the order of a|‘^^ for all t, 
as discussed in Sec. IIII A1 This makes the average over 
the trajectory time {x[^) to be w-dependent and differ¬ 
ent from the average over the Lyapunov time 
Applying the analysis without the division in regimes 


Sm leads to strongly fluctuating results (Figs. [nKa,c,e)). 
Much smoother results (Figs.EKb,d,f)) are obtained when 
we apply our method and compute only for t in 

the fully chaotic regime Sn- Looking at these smoother 
results we observe that the difference in Lyapunovs scales 
as 1/1 In ^1, but that even for w —>■ oo the sticky motion 
leads to a deviation from Daido’s relation 1®) (curves are 
shifted vertically). 


V. CONCLUSIONS 

In summary, we have proposed a method to charac¬ 
terize the dynamics of Hamiltonian systems with mixed 
phase-space based on time series of finite-time Lyapunov 
exponents. Using this method it is possible to define and 
study with high accuracy the time evolution of regimes of 
ordered, semi-ordered, and totally chaotic motion. This 
allows for an individualized characterization of the differ¬ 
ent stickiness mechanisms, improving alternative meth¬ 
ods based on the statistics of recurrence times or on the 
distribution of finite-time Lyapunov exponents. 

We applied our method to a chain of coupled stan¬ 
dard maps and showed how the frequency of different 
regimes - and the transition probabilities between them 
- are related to the volume of different phase-space re¬ 
gions. Using the consecutive time in distinct regimes 
we have reproduced previous results obtained using re¬ 
currence times and showed that our method allows for 
a significant improvement in the characterization of the 
sticky motion {e.g., in the determination of the scaling 
exponents). This indicates that our method can be used 
to characterize stickiness in general high-dimensional sys¬ 
tems and is particularly suited for cases in which different 
regions of sticky motion coexist. We have also shown that 
the dependence on the coupling strength of the largest 
Lyapunov exponents, after conveniently using our proce¬ 
dure, tend to follow only the qualitative universal prop¬ 
erties of fully chaotic system. 

Results obtained in a simple chain of standard maps 
confirm that our methodology can be applied to high¬ 
dimensional systems and problems of current interest, 
such as controlling Fermi acceleration [s^ , galactic mod¬ 
els [13, and plasma physics [13 ■ Another example of 
application is to associate each regime Sm with effective 
Hamiltonian functions, a procedure used to reproduce 
the complicated dynamics of kicking electrons |36| or the 
high harmonic generation in laser-assisted collisions [13 ■ 
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